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Abstract 

Quark mass corrections to the spin partonic structure function gi{x,Q^) and 
function F^{x,Q'^) are obtained at the order 0{as) along with the coefficient func- 
tions C^^^ and C^^) related to the Bjorken and Gross-Llewellyn-Smith sum rules. In 
the massless limit the difference between F3 and gi is encountered to be (as/7r)CF(l— 
x). The results for the functions C^"^^ and C*^^^ at m = agree with the previous 
MS-scheme calculations C(^) = C(^) = 1 - (3as/47r)CF. 
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Recently many theoretical efforts were spent to study radiative corrections to the 
polarized Bjorken sum rule (BSR) in deep inelastic eA^-scattering 
1 

J {gr{x,Q')-gr{x,Q'))dx = ^ 



where gA and gy are the constants in nuetron weak decay and C^^^ is the coefficient 
function of the axial current in the operator product expansion of two vector currents. 
Corrections to C^"^-* are known to coincide with those to another coefficient function C^^-* 
which is relevant for the Gross-Llewellyn- Smith (GLS) sum rule in neutrino-nucleon 
scattering 

1 

/ (F|^(x,Q2) + F3^^(a;,Q2))d^ = 3C(^). (2) 



The leading corrections to BSR and GLS were first computed in Refs. |0, ^ and 
appeared to be C^^^ = C^^^ = 1 — Ci?(3Q;s/47r). Next to leading order results for both 
sum rules can be found in Ref. [Q. Discrepancies between C^^^ and C^^-* arise only at 
order 0{a^) Q where 'light-by-light' disgrams appear. These calculations were performed 
in MS scheme with quark mass m = 0. From the other hand recent time mass dependent 
RG equations |^ and power corrections attract a great attention. Thus it is of interest 
to get mass dependence of coefficient functions. 

In the recent work there were computed 0{as) corrections to C^"^^ in on-shell 
scheme. It was noticed that there are contributions to the partonic structure function gi 
that survive when m 0. That is the diagram of the box type which is responsible for 
that. After momentum integration it develops an additional compensating factor 
which cancels with the mass in numerator resulting in finite terms. Thus the result for gi 
is different in dependence on whether mass m = from the very beginning or it is kept 
till the end of calculations. 

Below the partonic structure functions both gi and at order 0{as) with m ^ will 
be presented. Along this paper we use the on-shell renormalization scheme. 

First we note that the contribution of virtual gluon can be expressed in terms of 
renormalized elastic formfactors of currents. Corrections to the vector (axial) current are 
usually written as {q = p' — p) 

V, = u{p') {7,(1 + ^C^-^r(g)) + l-[q^^^]_^CpTUq)]n{p), (3) 

A, = u{p') |7,75(l + '^Cp:F^{q)) + l-{p + p')^^,'^Cp:F^{q)Y{p). (4) 
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These formfactors were computed earlier in QED and electroweak theory (see e.g. [|T0 
and references therein). We take them in the following form 

= - + 13^ - " 4(1-^^) ^ 

1 + ^ V 1 , 2.1 



•^2'' = -Y^log^' (6) 
•^i^ = -^r + Y3^1og^, (7) 



with 



^ v^rT4f- 1 

C7 : 



2 

m 



In the above formulae /i is a small 'gluon mass' (/i <^ t, m^) being infrared regulator and 
m is a quark mass. Formfactor JF^ is irrelevant and hence is omitted. The vector current 
is normalized such that J-'Y {q = 0) = 1. If such the condition is imposed then for the 
axial current one gets ^^(0) = 1 — (a;s/2vr). Formfactor JF^(O) = 1/2 is the well known 
anomalous magnetic moment. Let us emphasize that it is impossible to normalize both 
currents to unity in the presence of a nonvanishing mass. We shall return to this question 
later on. 

Using definitions (^,(§) the virtual contributions can be cast into the form 

1 a 



9i'— = 7;r^FS{l - x) {Sj^r + 4-^2^} , (9) 



-.virtual 

~ 2 471 

^virtual ^ i^C^5(l - a;) {atY + 4-^i^} ■ (10) 

Next we turn to the real gluon contributions. Again one faces with IR divergencies 
which are due to soft or collinear gluon emission. As above they are regularized by letting 
a gluon have a small nonzero mass /x. Calculating strightforwardly the diagrams one 
obtains functions which have a nontrivial dependence on /i. It causes singularities like 
1/(1 — a;) at the end point of x-integration (in presence of a small mass fi x is varying 
in the interval < x < 1 — 2r/i/m). Of course after integration over variable x the 
fi dependence is cast into that as in formulae (P|),([T(]|) with the opposite sign so that /i 
cancels in the whole moments. To cancel out fi before integration we rewrite result using 
the well known 'plus distribution' |]11[, i.e. every function /(x) being singular at x = 1 is 
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replaced by 



1— 2r/i/m 

f{x) = U{x)+6{l-x) I f{z)dz, 
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where the hmit — > is imphed. After some transformations the structure functions can 
be presented as 



1 as 



^ (a;) = t^itCf 5(1 - x)Tl + 4(1 + 2r) 
2 47r 



L 



1 — X 



;i-x). 



log 



1 + 2rx + Vl + 4rx2 



+ C1 + C2L , (12) 



(13) 



a/1 + 4rx^ 1 + 2ra; — Vl + 4rx^ 
Here iy(x) stands for either gi or F3. All infrared divergencies now are absorbed in 
coefficients 71. By explicit calculation we have found for 71, Ci and C2 



7^9l = ^SJ^^ - 4J^2 - 4 + 8 logr - 2 ^ 



5^2 + 4^ + 5 



02 



log^ 



1 + 02 



T^Fa = _4j-v _ _ 4 ^ 8 log^ _ 10—^ log e 

1 — 0^ 

^ ^[l~x){3 + 6x-8x'^) 



(14) 
(15) 
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(1 + 4rx^)2(l — X + rx)2 
+ 4rx(l - a;)(2 + 15x - 15x^ + 4x^) 

+ 4r2a;2(l + 24x - 29x^ + 6a;^) + 8r V(5 - a;) 



1 



/3 ^ 

(1 + 4ra;^) (1 — x + ra;)^ 
+ 2rx{7 - 3a; - 4a;^) + 8r^x^ 



1 -x)(7-9a;-2a;-4a;^ 
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1+X + 2r(2 + x + 9x^ - Ax^) 



(l + 4rx2)2 
+ 4r^x^(8 + 5x - x^) + 64r^x'^ 



^2 



1 + X + 4r(l + x) + 16r^x^) 



(16) 

(17) 

(18) 
(19) 



1 + 4rx2 

From the formulae (P),(p!0|) and (p^-(p^) one can see that logyu's cancel as well as Ti 2^ 
and we are left with IR safe expressions. 

If ^ then the formulae for gi, F3 greatly simplify. Adding the Born contribution 
and keeping only leading terms in w? jQ^ we arrive at 



1 + X^ 



[\-x\ 



'log(x(l — x)) 

1 — X 



^ +2 + 8x + 2(l + x)log(x(l-x))} (20) 
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F3(.) ^ 1^(1 - X) + 1^C.{2 [Isil - .) + ^] log I, - 5.(1 - X 



-4 ( log(^(l ^ ^)) j +6 + 4x + 2(l + x)log(x(l-x))} (21) 

(22) 



(1 — x)_|_ Y 1 — X 



As it was mentioned in the beginning of the paper corrections to the structure function 
F3 evaluated in massless theory coincide with those to gi. Here we see from (pOD,(pTD that 
this is not true in massless limit of the massive formulae. The difference is totally due to 
coefficients ci's. It is worth noting that the result obtained should not depend on an in- 
frared regularization procedure provided quark mass m tends to zero after IR singularities 
are canceled out. Eqns. (^),(p]) yield that gi develops an extraterm — Ci?(a;s/27r)(l — x). 
In fact this term defines a fermion helicity-flip probability P+_ (x) and was studied in Ref . 

0- 

Using formulae (|I^)-([l9D we obtain first moments of gi and F3 with no approximation 
made 



1 



/Ql 
wix) dx = 1 + —Cf5, wix) = gi(x), F^x). (23) 
47r 



Where the factors 5's read 

. l-16r^, , 2-17r + 16r2 i ^ - 28r\ ^ 

S'' = -4 + -^^Hr) — ^ logr + log^, 24 

2 2r(l — r) rV 1 + 4r 

5-4r ^ ^, , l + 6r- 16r2 + 8r3 1 + lOr + 24r2 

d ^ = 8rI[r)-\ logr ^ — logfe*. (25) 

1 — r r(l — r) rvl + 4r 

Function /(r) can be expressed through Euler dilogarithm function Li2 and has the fol- 
lowing representation 



^ , , , dx , 1 + 2rx + vT+4rx^ 

I[rj= I = log . : 

vl + 4rx^ 1 + 2rx — v 1 -|- 4rx^ 

^ :{Li2(-t) - Li2(t) + Li2(-at) + 112(7) - ^Li2(a2) + ^ + l log' t}, (26) 







2y/¥' ' ' ' ' 't' 2 ' ' 4 2 

a= (1- v^)/(l + v^), t = VTT^-2^. (27) 

In high region it turns to be /(r) = 2 — logr + r(2/9 + (5/3) logr)+0(r^ logr). 

Eqns. (p3|)-(^) together with (|26|) give the values of the first moments of the partonic 
structure functions in the whole region of > 0. In deep inelastic case [m? /Q"^ —>■ 0) we 
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obtain gf^ = 1 — Cir(5«s/4vr) while within the massless approach it is 1 — Ci? (3as/4'7r). 
This resuh was found earher in Ref. 

Let us consider now the coefficient functions. The OPE says that 

2g?{Q')=C^^\Q^)A^\ (28) 
2FP{Q^)=C^^\Q^)A'^^\ (29) 

where As are defined from operator matrix elements 

(p,s|V^(g)757^V^(g)|p',s) =M(p)7^75M(p')^^^^(t), (30) 

{p\i'{q)lM<i)\p')=^ip)lAp')A'^'\i)^ (31) 
t = {p-p)\ (32) 

with quark field if). Matrix elements (|30| ) ,(|3T| ) must be taken at zero momentum transfer 
t = in the on-shell scheme. For the vector current the identity A^^^ = 1 + (a;s/7r)jF|"(t = 
0) follows and it is equal to unity because of the current normalization while for the axial 
we have A^^\t) = 1 + {a^/T:){J^^(t) — J-'Y {t)). The structure in the parentheses looks like 

^i^W--^rW = Y^log^, (33) 

with 9 defined as in (|^) and r = rn^/t. For large t (|33| ) vanishes as it should be due to 
chiral invariance. This situation is realized when m is identically equal to zero, when the 
limit t ^ corresponding to the forward matrix element could be easily taken. However, 
if one take this limit before setting m = the mass terms come into the game and (p^) 
becomes —1/2. One can check this using (^-(H). As a result C^"^-* differs form the g'l'' at 
m = by a finite term 

= (l - ^C'.) (l - |C.)" = 1 - I^C. + 0(^). (34) 

Let us summarize now the results. Coefficient functions C's with a nonvanishing 
fermion mass look like 

C^A) ^ ^ ^ ^Cf{S'' + 2), (35) 

4:71 

C(^) = 1 + ^Cf6^^ (36) 

4:71 



where functions 5 are given by (plD , (psf) . In Fig.l there are the plots of the corrections to 
coefficient functions C-^^ C^^^ versus m? /Q"^. In the deep inelastic limit both corrections 
coincide with each other in agreement with the values quoted in literature. 
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Up to the terms 0{m? /Q'^Y the Eqns. (|3 



1 



a 



read 



10\ 

3 + 7)^(t^°Sq^-y) 



Air 



3 + 



(slog 



+ 4 



(37) 
(38) 



The discrepancy between C'-^^ and C*^^^ manifests a violation of the Crewther relation 
[p^ by mass corrections. 

The magnitude of the strange quark mass correction is about 0.12 of the massless 
one- loop result at = 2 GeV^. For the light quarks, the mass contribution is negligible 
if one uses the current quark mass of order of few MeV. If, however, one takes into account 
that such a scale should be non-observable and substitute instead the scale of order of a 



pion mass |14], the result is still about 0.1. 



The calculated corrections are in fact the first example of the NLO mass dependence 
in QCD. Generally speaking this would require to calculate the 2-loop mass- dependent 
anomalous dimension and one-loop coefficient function. However in the case at hand the 
anomalous dimension is zero and the calculated contribution provides the final result. 



After this work was completed, the paper |]15[ appeared, where the coefficient functions 
for heavy quarks are investigated in the limit > m?. 

We are indebted to A.L. Kataev, S.V. Mikhailov and I.V. Musatov for useful discus- 
sions. O.T. is grateful to J. Collins for elucidating correspondence and to W. van Neerven 
for valuable comments. 
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Fig. 1. Corrections to the Bjorken (solid line) and Gross-Llewellyn-Smith 
(dashed line) sum rules versus vr? jQ'^ . Coefficient functions are written in the 
form = 1 + (as/47r)CFA^'^. 
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